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All-oxide magnetic tunnel junctions (MTJs) incorporating functional materials as insulating bar- 
riers have the potential of becoming the founding technology for novel multi-functional devices. We 
investigate, by first-principles density functional theory, the bias-dependent transport properties of 
an all-oxide SrRuOa/BaTiOa/SrRuOs MTJ. This incorporates a BaTiOs barrier which can be found 
either in a non-ferroic or in a ferroelectric state. In such an MTJ not only can the tunneling magne- 
toresistance reach enormous values, but also, for certain voltages, its sign can be changed by altering 
the barrier electric state. These findings pave the way for a new generation of electrically-controlled 
magnetic sensors. 



The control of the spin dependent tunneling between 
two ferromagnetic electrodes separated by an insulating 
barrier has enabled enormous advances in many mag- 
netic data storage technologies, in particular since ex- 
tremely large tunneling magnetoresistance (TMR) was 
measured. The progress in producing magnetic tunnel 
junctions (MTJs) with large TMR was initially limited 
by the use of amorphous tunnel barriers. The situation 
however changed after the prediction [l], 0] and subse- 
quently experimental realization 0,0] of epitaxial MTJs. 
Since then, room temperature TMR in excess of 600% 
has been demonstrated in MgO-based devices [Bj]. 

In general, for amorphous barriers the spin polarization 
of the tunneling current and hence the TMR magnitude, 
depend solely on the electrodes' density of states (DOS) 
at the Fermi level, Ep [6j. In contrast, perfectly crys- 
talline tunnel barriers are wave-function symmetry se- 
lective and make the tunneling process sensitive to their 
electronic structure. As a result the amplitude and even 
sign of the TMR may depend on the barrier itself. The 
understanding of such a concept suggests that one can 
engineer the TMR by carefully selecting the insulating 
barriers to be epitaxially grown on magnetic electrodes. 
Ferromagnets Q and ferroelectrics [1, |J are of particular 
interest as functional barriers. 

Ferroelectric materials possess a spontaneous electric 
polarization whose direction can be switched by an elec- 
tric field. This makes ferroelectric-based MTJs fully 
multi-functional devices able to respond to both electrical 
and magnetic stimuli d, [l(J ■ Importantly ferroelectrics 
can be grown epitaxially on a variety of substrates 11 1 
but in particular on other oxides. Since epitaxial growth 
is a pre-requisite for large TMR, the prospect of all-oxide 
junctions appears particularly attractive. Such a type of 
MTJ is investigated in this Letter. We demonstrate the- 
oretically a huge TMR and more importantly we show 
that the TMR sign can be reversed with bias, at a criti- 
cal bias which depends on the ferroic state of the barrier. 
Our results are rationalized in terms the band-structure 
match between the ferroelectric insulator and the ferro- 



magnetic electrodes. 

Density functional theory (DFT) calculations are per- 
formed with the local basis set code siesta Struc- 
tural relaxation is obtained with the generalized gradient 
approximation (GGA) of the exchange and correlation 
functional [l3| • This gives a satisfactory device geometry, 
but it produces a rather shallow band alignment mainly 
because of the DFT-GGA gap problem. In order to make 
up for this shortfall the electronic structure used for the 
transport calculations is that obtained with the atomic 
self-interaction correction (ASIC) scheme 14[, which im- 
proves drastically the electronic properties of both bulk 
BaTi03 [HI and SrRu03 [3] . Unfortunately the approx- 
imate ASIC energy functional is not sufficient to produce 
good structural parameters and in particular the BaTiOs 
ferroelectric state cannot be stabilized. This is a current 
limitation of the method, which otherwise has been suc- 
cessful in predicting the electronic properties of oxides 
18]. For this reason we perform ASIC transport calcu- 
lations at the GGA relaxed structural parameters. For 
all the calculations we use a 6x6x1 fc-point Monkhorst- 
Pack mesh and a grid spacing equivalent to a plane- wave 
cutoff of 800 eV. 

Electron transport is computed with the SMEAGOL 
code (l9l [2(j |. which combines the non-equilibrium 
Green's function scheme with DFT. Since SMEAGOL in- 
terfaces siesta as the DFT platform, we employ here the 
same parameters used for the total energies calculations. 
In brief, the total electronic current is given by 



(1) 



where a labels the spin (t,|), T a (E; V) is the energy de- 
pendent transmission coefficient for the bias V, Jl/r is 
the Fermi distribution function evaluated at E — ^l/r 
and fiL/R = Ep ± ^- is the chemical potential of the 
left/right electrode. If the junction is perfectly transla- 
tional invariant in the plane orthogonal to the transport 
direction, T a is obtained by integrating the fc-dependent 
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TS over the 2D Brillouin zone of volume fiez , 



T a (E,V) 



BZ JBZ 



dk T2(E; V) 
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We initially perform relaxation of bulk BaTiC>3 and 
SrRu03 under an in plane compressive strain emulat- 
ing the common epitaxial growth on SrTi03. The re- 
laxed cells are then used to construct the transport super- 
cell, which comprises of six BaTiOa unit cells (~2.5 nm) 
sandwiched at either side by three SrRu03 ones. The 
SrRu0 3 /BaTi0 3 interface is SrO/Ti0 2 , due to the ex- 
perimentally observed volatility of the R11O2 termination 
[2l| . We consider two structures. In the first non-ferroic 
(NFE) structure the atoms are frozen artifically in their 
ccntro-symmctric positions with the interfacial distance 
given by an average between the BaTiOs and SrRuC>3 
c-lattice constant. In the second the supercell is further 
relaxed with respect to the atomic coordinates, resulting 
in a stable ferroelectric ground state (FE structure). 
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d xz and d yz ). Finally the d x 2_ y 2 and d xy states have A2 
and A 2 < symmetry respectively. Importantly an incident 
Bloch state in the electrodes can couple to a given evanes- 
cent state in the insulator, and then sustain a tunneling 
current, only if the two share the same symmetry. 

The left panel of Fig. [5] shows the SrRuOs band- 
structure close to Ep along the direction of the trans- 
port. At Ep only a doubly-degenerate minority A5 state 
is available, in contrast to previous DFT calculations, 
where both minority Q,) A5 and majority (j) Ai bands 
were found [lj|. Such a discrepancy is due to the use of 
the GGA functional in Ref. J10|, which underestimates 
the Ru d manifold exchange splitting [lj|. Note that a 
large spin splitting is expected based on point contact 
Andreev reflection experiments 
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FIG. 1: (Color online) Relative atomic displacement for 
Sr, Ba (open circle) and Ru, Ti (filled circle) in the 
SrRu0 3 /BaTi0 3 /SrRu0 3 MTJ investigated. The displace- 
ments are with respect of the O atoms in the same plane, 
5 = (^cation — zo), with z being stack direction. The device 
geometry is presented in the lower part of the figure. 

The atomic relaxed displacements, <5, with respect to 
the planar O positions are shown in Fig. [TJ At the cen- 
ter of the the BaTiOs slab Ti displaces by 0.14 A, which 
is significantly smaller than the value of 0.23 A of bulk 
BaTi03 experiencing the same strain. Note that GGA 
overestimates the volume and atomic distortions associ- 
ated with ferroelectricity in BaTiOs resulting in a "super- 
tetragonal" structure. Such an overestimation, while re- 
sulting in a polarization greater than the experimental 
one, has not a significant qualitative effect on our re- 
sults. The interfacial SrRu03 layers, as expected, also 
contribute to the polarization [221 ] . 

The symmetry of the electronic bands of both the fer- 
romagnetic electrodes and the insulating spacer dictates 
the transport properties. A wavefunction, whether prop- 
agating or evanescent, is described in terms of irreducible 
representations of the crystal's symmetry group. For a 
cubic space group, the Ai symmetry transforms as a lin- 
ear combination of 1, z and 2z 2 — x 2 — y 2 functions, while 
the A5 as a linear combination of zx and zy (e.g. p x , p y , 



FIG. 2: (Color online) ASIC-calculated band-structure along 
the transport direction (F — > Z) for centro-symmetric tetrag- 
onal SrRu0 3 (left), and BaTi0 3 (middle) both in the NFE 
(green) and FE (red) configuration. The wavefunction sym- 
metries of the bands close to Ef are indicated. The right 
panel reports the complex band-structure for BaTi0 3 . The 
energies are aligned with the Ef of SrRu0 3 . 

In the right panel of Fig. [5] we plot the BaTiOs real 
and complex band-structure. In contrast to MgO, where 
states with A5 symmetry decay significantly faster than 
those with Ai Q|, in NFE BaTi0 3 the Ai and A 5 sym- 
metries have comparable decay rates. In particular close 
to the valence band top the slower decay rate is for Ai, 
while the situation is reversed at the conduction band 
minimum. The enlargement of the bandgap associated 
with the FE order results in an increased decay rate for 
all the symmetries. The effect is more pronounced for 
A 5 close to the top of the valence band where now the 
Ai symmetry primarily contributes to the tunnel con- 
ductance. 

We begin our analysis of the transport properties from 
the NFE structure by showing T(E) at zero bias for the 
parallel (PA) and antiparallel (AP) magnetic alignment 
of the electrodes (Fig. [3]). In the PA configuration T(E) 
close to E-p is dominated by the minority spin channel. 
This is expected from the band-structure of SrRu03, 
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which presents only a A5 4- band along the transport 
direction for energies comprised between —0.8 eV and 
+0.1 eV. The minority conductance in this energy range 
is five orders of magnitude larger (T^ ~ 10 -7 ) than that 
for the majority spins. For E > 0.1 eV there is a sharp 
rise in T^, due to the Ai f band now contributing to the 
conductance. In the energy window 0.3 eV < E < 0.8 eV 
there are no minority states available and drasti- 
cally drops. A similar drop, due to the lack of minority 
SrRuOs bands is found at -2.5 eV. In the AP configura- 
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FIG. 3: (Color online) Transmission coefficients, T(E), as a 
function of energy for the NFE structure. The middle panel 
is for the parallel magnetic configuration, while the lower one 
is for the antiparallel. At the top we report again the SrRu03 
band-structure at the T point of the 2D transverse BZ. The 
dotted line at eV denotes Ef- 

tion the electron transmission occurs between majority 
(minority) states in the left-hand side electrode and mi- 
nority (majority) in the right-hand side one, so that T(E) 
for both the spins (identical) is essentially a convolution 
of those for the majority and minority spin channels of 
the PA state. As a consequence there is a drastic sup- 
pression of T(E) in the regions -0.8 eV < E < 0.1 eV and 
0.3 eV < E < 0.8 eV, where respectively the A 5 I and 
Ai f bands in one electrode are not paired in the other. 
In particular T(E-p) for the AP configuration is orders of 
magnitude smaller than in the PA one. Note that our dis- 
cussion is based on the band-structure at the T-point of 
the 2D transverse BZ, for which the decay is the smallest 
and the transmission the largest. However, also Bloch 
states with larger transverse wave-vector contribute to 
the transport and produce a residual transmission. 

The spin-polarized current for both the PA and AP 
configurations and for both the NFE (top panel) and 
FE (middle panel) structures are shown in Fig. |U where 
we focus on the low voltage region in which the cur- 
rent is due entirely to tunneling (the broader I-V are 
displayed in the insets). The most distinctive feature 
emerging from the I-V curves is the presence of nega- 
tive differential resistances (NDR) for the PA alignment, 
originating from the movement of the Ai j" band-edge 
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FIG. 4: (Color online) Total current per unit area, I, as func- 
tion of voltage, V, for the NFE (top panel) and FE (middle 
panel) structures. In the bottom panel we present the TMR 
as a function of voltage for both the geometries. In the insets 
the I-V is presented over a larger current range (the units are 
the same as in the main figure). Note that at the on-set of the 
BaTiOa conduction and valence bands the current increases 
by three orders of magnitude over its low bias value. 



with V. Because of the NDR the relative magnitude of 
the current for the parallel (/ PA ) and antiparallel align- 
ment (I AP ) can be reversed, i.e. the TMR changes 
sign with V. This is demonstrated in the lower panel 
of Fig. |4l where we present the "pessimistic" TMR ra- 
tio, TMR=(/ PA - J AP )/(/ PA + 7 AP ), as a function of 
bias. Clearly TMR sign inversion is observed for both 
the NFE and the FE junctions for voltages in the range 
0.7-0.9 V. Furthermore for V ~ 0.7 V the TMR for the 
NFE junction is positive, while that of the FE one is 
negative, meaning that subtle changes in the barrier elec- 
tronic structure, such as those induced by ferroelectric- 
ity, are sufficient to change the sign of the TMR. Note 
also that the TMR values reported here are actually ex- 
tremely large. For instance for both the NFE and FE 
junctions and voltages \V\ < 0.4 V the optimistic TMR 
[(/ PA - / Ap )// Ap ] is around 5,000%. 

The I-V curve can be rationalized by looking at the 
dependence of T(E) on the bias [23J, which is presented 
in Fig. [5] This is mainly determined by the shift of the 
electrodes' Ai and A5 band-edges with V: for positive 
voltage the band-structure of the left electrode is shifted 
by +eV/2 (e is the electron charge) and that of the right 
one by — eV/2. At a given energy a large T is found 
only if a band of the same symmetry and spin is found 
in both electrodes at that energy. For PA alignment at 
V = 0, the minority spins dominate the transmission up 
to 0.3 eV, after which one encounters the A^ upper band- 
edge and is drastically reduced. As V is applied, the 
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Ag band edge is shifted to lower energies in the right 
electrode (for V > 0), so that for V — 0.6 V the high 
transmission region extends only up to Ep, and for V — 
1.2 V it extends only up to Ep — 0.3 eV. This is the origin 
of the NDR found for the PA alignment. In contrast 
for the AP configuration T is small for energies below 
0.1 eV after which it drastically increases because of the 
Ai conduction bands (see Fig. [3]). With increasing V the 
Ai band in the right electrode is shifted to lower energies, 
so that there is a rather large transmission inside the 
bias window, and eventually the AP current therefore 
becomes larger than the PA one. This results in the 
TMR sign change at about 0.7-0.9 V. 




1 1 1 1 1 1 1 
i i i i i i i 


111111/ 


■ 




■ Is' 






1 1 i h„i i'V 


1 1 1 | 1 r 1 






., Hi , , , 





-1 -0.5 0.5 
£(eV) 



1 -1 -0.5 0.5 
£(eV) 



FIG. 5: (Color online) Transmission coefficient, T(E;V), as 
a function of energy and for different bias voltages. The left 
panels are for the parallel configuration and the right ones for 
the antiparallel. In the same figure we report data for both 
the NFE structure (t black solid line, 4- black dashed line) 
and for the PE one ft red solid line, J, red dashed line). The 
vertical dashed lines mark the borders of the bias window. 

The main effect of the ferroelectric order on the trans- 
port is an increase of the BaTiOa band-gap, i.e. an in- 
crease of the Ai and A5 decay coefficients (see Fig. [2]). In 
particular, states with A5 symmetry decay significantly 
faster in the FE MTJ with respect to the NFE one. This 
results in a global reduction of the transmission although 
other general features remain rather similar in the two 
cases. The comparison between T(E; V) for the FE and 
NFE junctions is also presented in Fig. [SJ Below Ep 
one may note a substantial reduction of the transmission 
when going from NFE to FE for both PA and AP align- 
ment as a consequence of the increased A5 decay rate. 

In conclusion, we have demonstrated huge TMR in an 
all-oxide ferroelectric MTJ, the sign of which can be in- 



verted as the applied bias increases. Furthermore the sign 
inversion occurs at different voltages for different ferroic 
states of the barrier. Our finite-bias results are explained 
in term of the electrodes and the barrier band-structures. 
The possibility to control the TMR by manipulating the 
ferroic state of the barrier in an MTJ opens a potential 
avenue for the electrical control of magnetic devices. 
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